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THERMODYNAMIC BEHAVIOR OF Ga 3+ , Al 3+ AND Fe 3 + 
DISTRIBUTIONS IN GARNETS 

ABSTRACT 

Data on the distribution ol Ga 3+ , Al 3+ and Fe 3+ among the tetrahedral 
(d) and octahedral (a) sites in certain garnets were examined in terms 
of a thermodynamic model which yields the following expression for the 
distribution constant K°: 


R , X d - V 

X a (1 * X cP <*ad 

Here the X's are site occupancy fractions and the f's are equivalent to 
activity coefficient products. It is shown that for the elements discussed 
(i^/tad) ~ N corresponding to ideal mixing on the individual sites. Thus 
we obtain 
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Ga 3+ 

(Ga 3+ + 

Fc* 3 + ) 

6.71 

1643 

Y Ga Al. 0,, 

3 x 5 x 1 2 

Ga 3+ 

(Ga 3+ + 

Al 3+ ) 

-2.24 

- 

Eu, Ga Ft. 0 

3 x 5 x 1 2 

Ga 3+ 

(Ga 3+ + 

Fe 3+ ) 

4.46 
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where X ranges from zero to five. 


v 



With the presently available precision the calculated temperature 
sensitivity of K° is such that variations ol' several hundred degrees in the 
sintering temperature should be apparent 

It is stressed that this form ol ideal mixing on the individual sites 
corresponds to nonideal mixing of the component end members V 5 F«? 5 
0. 2 and Y 3 Ga s Oj , when referred to their mole fractions and mean molar 
volumes. This is consistent with the apparent positive deviation of the 
mean molar volume from a straight line. 


vi 



THERMODYNAMIC BEHAVIOR OF Ga 3 + ,AI 3T AND Fe 3+ 
DISTRIBUTIONS IN GARNETS 


INTRODUCTION 

Nearly all crystalline materials are of a nonstoichiometric character and at 
the same time have atomic structures in which each atom or ion is distributed 
over two or more nonequivalent structural positions. In the garnets Ga 3+ , Al 3t 
and Fe 3+ are distributed among two different crystallographic sites with octa- 
hedral (a) and tetrahedral ( d ) coordinations of the O ' 2 ions. The nature of the 
distributions or oraer-disorder relations between these positions are functions of 
the total crystalline composition, the temperature and the pressure in the 
equilibrium state. In a disequilibrium state the distributions are also dependent 
on the previous history of the sample and on the time and provide an interesting 
application of chemical kinetics. 

The first comprehensive treatment of the thermodynamic theory for order- 
disorder phenomena was provided by Bragg and Williams 1 for the simple binary 
alloys. Analogous treatments were subsequently given by others for the spinels. 
Especially noteworthy contributions were made by Neel , 2 Smart 3 and by Callen, 
Harrison and Kriessman . 4 Recently Borghese 5 interpreted the magnetic proper- 
ties of certain spinels and garnets in terms of the model in which mixing on the 
individual lattice sites has the character of ideal solutions. Analogous equilibrium 
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equations were also derived for the distribution of Mg 2+ and Fe 2+ among the 
nonequivalent structural sites of pyroxenes and amphiboles. 6 As distinguished 
from the previous treatments of metallic alloys, spinels and garnets, the latter 
equations were formulated in terms of the atomic fractions (site occupancy frac- 
tions) of the cations on the individual sites, a procedure which will be adopted 
in this paper also. This formulation has the advantage of yielding more simpli- 
fied data plots. 

Although the effect of temperature on the int.racrystalline distributions forms 
an important part of the problem and is implicit in all the theoretical treatments, 
this aspect of the problem has not been sufficiently appreciated in much of the 
experimental work on garnets so that irmre often than not the temperature of 
preparation of the specimens is not given. We shall attempt here to analyse 
the character of the solid solutions and to calculate the effect of temperature on 
the distribution isotherms. 

OBSERVED DISTRIBUTIONS 

There are three sets of experimental data for the distribution of Ga 3 * and 
Fe 3+ among the tetrahedral and octahedral sites of the garnet Y ? Ga x Fe s _ x ,, 
where x can in theory range from zero to five, corresponding to complete sub- 
stitution of Ga 3 * for Fe 3 ’. These data were obtained by the analysis of residual 
magnetic moment, 7 nuclear magnetic resonance 8 and by X-ray and neutron 
diffraction. 9 In Fig. 1 they are plotted in terms of the atomic fraction Ga 3+ / 
(Ga 3 * + Fe 3 *) (where the chemical symbols represent numbers of moles) of 
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gallium on the tetrahedral and octahedral sites of the total Ga 34 + Fe 3+ on 
each site. 

Data for the distribution of Al 3+ and Ga 3+ in Y, A1 Ga r O, „ available 

3 x 5 -x 12 

from X-ray diffraction experiments 10 are given in Table 1. Magnetic data de- 
rived from the curve of Geller et al 11 for Y, Al x Fe s x O, , and as discussed 
by Borghese 5 are shown as points in Fig. 2 for comparison. 

Table 1. Distribution of Ga 3 * and A 3 * in Y, Ga Al , O . _ 
and the corresponding equilibrium constant of exchange. 

Data are taken from Marezio et al. 10 


x 


Ga 3+ /(Ga 34 

+ Al 34 ) 

K° 

Tetrahedral (<1) 

Octahedral (n ) 

0.290 

0.080 

4.09 

0.491 

0.223 

3.30 

0.053 

0.470 

2.12 

0.803 

0.744 

2.17 


Some experimental points obtained for the distribution of Ga 3 + and Fe 3+ in 
the garnet Eu 3 Gn x Fe g x O j , by Mbssbauer spectroscopy are shown in Fig. 3. 

Unfortunately it is not possible to evaluate some major uncertainties in the 
experimental site occupancy fractions because of lack of information on the con- 
ditions of sample preparation. However according to Fischer et al 9 the X-ray 
and neutron diffraction data embody probable errors of about 4 to 5 percent of 
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the measured values. While the scatter of points in certain parts of Fig. 1 is 


substantial, the overall agreement between the three sets of data is quite good. 

VARIATION OF THE MEAN MOI \R VOLUME WITH COMPOSITION 

The variation of the unit cell volume of Y Ga Fe c O, with Ga 3+ /(Ga 3+ 

3 x 5 -x 12 ' 

Fe 3+ ) is shown in Fig. 4. These data, which have a direct bearing on the thermo 
dynamic interpretation of the distribution data are taken from the work of Gilleo 
and Geller, 13 Espinosa, 14 Euler and Brace 15 and Fischer et al. 9 The chief 
interest here is in the possible deviation from a straight line between the molar 
volumes of Y 3 Ga s 0 J2 and Y 3 Fe g 0 J2 since adherence of the data to such a 
line would be required for an ideal solid solution between these macroscopic 
components. Unfortunately the uncertainties in the cell volume data do not allow 
this question to be answered unambiguously. However, it is clear that even if a 
straight line were drawn through the highest volume point (that of Euler and 
Bruce) for Y 3 Ga s 0 J2 , most of the intermediate points would still lie above this 
line. From the most probable straight line (drawn as dashed) this result is even 
clearer. If this is the true picture, is is clear that there is a positive excess 
molar volume and that the garnet solid solutions as a whole is non -ideal in 
character. 

THERMODYNAMIC CHARACTER OF THE SOLID SOLUTIONS 

If the dependence on composition of the free energy of exchange between the 
noncquivalent lattice sites is largely confined to the entropy of mixing terms, 
it is convenient to write 16 the following expression for a quasibinarv crystal 


such as a garnet: 
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( 1 ) 


K, n - X a ) / 


del 


X., (1 - X.) t 


ad 


Here the X 's are the site occupancy fractions of one of the two exchangable species 
on the octahedral *a) and tetrahedral (d) sites, K° is the equilibrium constant and 
ihe 0 factors arc* analogous to activity coefficient products and are derived from 
kinet*c formulae. By definition K° is a function only of the temperature and the 
pressure (or volume) while the 0's are also functions of the total composition. 
However, only the temperature dependency is significant for these condensed sys- 
tems at low to intermediate pressures. 

The site occupancy fractions are related to the total mole fraction \ t of a 
given end member (such as Y 3 Ga_. O j2 ) as 



The simplest and most convenient form of equation (1) is one for which ( t (i J 
4 <1(j ) = 1 and in which all compositional dependency is confined to the entropy 
of mixing on the individual a and d sites. This then is a form of "ideal" solu- 
tion, which is however sharply differentiated for the ideal solution lor the crystal 
as a whole and for which the activity n of a component end member (such as 
Y 3 Ga s 0, 2 ) is set equal to X t . Rather this activity (normalized to 1/5 Ga s 
0 12 ) is 



Obviously (2) and (3) can be equal only if X d = X a when K° = 1. However in the 
model under discussion K° ^ 1 even at the highest temperatures except as a 
limiting case since, unlike the nonequivalent sites in metallic alloys, the 
tetrahedral and octahedral sites always retain their identities. 

The curve of Fig. 1 is a fit of equation (1) to the observed data for Y 3 Ga x 
Fe 5 _ x 0 12 . In this case X d and X a refer to the atomic fraction Ga 3+ /(Ga 3+ + 
Fe 3+ ) on the tetrahedral and octahedral sites respectively. Except for some 
local scatter the data are quite well represented by the curve corresponding 
to K° = K ’ = 6.71. 

Fig. 2 shows the analogous curve as derived from equation (1) with X d and 
X n equal to the atomic fraction A1 3+ /(A1 3+ + Fe 3+ ) on the tetrahedral and octa- 
hedral sites of Y 3 Al^ Fe. x 0 ]2 . This curve corresponds to an equilibrium 
constant K° - K~ = 3 as derived by Borghese. 5 Borghese has shown that his 
version of equation (1) combined with his formula for the magnetic moment yields 
a better representation of the magnetic data than do the plotted points in Fig. 2. 

In fact it may be shown that the curve from which the latter are derived is incon- 
sistent with the thermodynamic model and equation (1) in any form. We therefore 
accept the value K i = 3 as derived by Borghese. 



The data of figures 1 and 2 should be directly related to those of Marezio 
et al as given in Table 1 if certain simplifying assumptions are introduced, for 
we have the following related exchange reactions and equilibrium constants: 

Ga (a) + Fe (cf) ^ Ga ( d ) + Fe (a) Kj 

Al (a) + Fe ( cl ) Z Al ( d ) + Fe (a) K° 

Ga (a) + Al (c/) 2 Al (a) + Ga (d) K' 

in which a and d refer to the octahedral and tetrahedral sites. If we now assume 
that the mixing properties of the ions are invariant with respect to which other 
ion is present we obtain 



K° = 2.24 


Table 1 shows the values of K° as derived from the data of Marezio et al . in 
Although the correspondence with K° = 2.24 is quite good for several specimens 
the others differ considerably from this value. It is not known if these differences 
are attributable to the uncertainties in the data such as arise from the conditions 
of sintering or to a failure of the assumptions employed in the model. 

The distribution curve and the corresponding equilibrium constant KA for 
the distribution of Ga 3+ and Fe 3+ among the tetrahedral and octahedral sites 
of Eu 3 Ga x Fe 5 _ x 0 12 are shown in Fig. 3. In this case the curve gives a quite 
good fit with K° = 4.46 if the most Ga-rich sample is disregarded. It is clear 
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from a comparison of these data with Fig. 1 Hiat the substitution of europium 
lor yittrium tends to favor Fe J + in the tetrahedral site. However here again it 
must be remembered that the effect of the thermal history is unknown. 

Taken together, figures 1 and 1 point up an interesting feature of the thermo- 
dynamics. If the data of Fig. 1 truely correspond to the simplified model with 
( ) - 1 the system is a form of ideal solution in which the only compositional 
dependence of the free energy comes from the entropy of mixing on the individual 
sites as given by equation (3). However this result in turn implies non-ideal 
mixing between the macroscopic components Y 3 Ga. O r and Y 3 Fe 0 J2 in which 
the sites are iot distinguished. The apparent positive deviation from a straight 
line of the me 'n molar volume in Fig. 4 is compatible with this interpretation. 

EFFECT OF TEMPERATURE 

Of primary interest is the determination of the effect of temperatures of 
annealing on the equilibrium constants. This effect could of course be evaluated 
directly if H° the standard entholpy of exchange were known, since 

1 n K’ H° 

T IF 

However, in the absence of independent knowledge ol II 0 we can estimate G°, 
the Gibbs free energy of exchange from the distribution data of Fig. 1 for which 
the annealing temperature is known, since 
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AG° = - RT 1 n K 


( 4 ) 


According to Fischer et al 9 the temperature of sintering for their specimens was 
1643°K and since = 6.71 we obtain A G° = - 6220 calories. Also we have the 
general relation 


3AG° 


where AS 0 is the standard entropy of exchange. 

Now we have reason to suspect that A S° will be small since the system is 
condensed and no large volume changes are involved. Also, this anticipation 
has recently been substantiated for the case of the order-disorder relations of 
Fe 2 ^and Mg 2 ^ in (Mg, Fe) Si0 3 (orthopyroxene), U ' which represents a situation 
closely parallel to the garnets. 

Thus if ASj 0 we have 

AHj - AG° . 

and the temperature effect can be estimated directly for equation (4). This has 
been done for five isotherms in Fig. 5. 

If Fig. 5 is valid indication of the temperature sensitivity of K° , it is clear 
that from the amount of scatter shown in Fig. 1 variations of several hundred 
degrees in sintering temperatures should be detectable. It is unfortunate that the 
thermal histories of most of the specimens are unknown. 
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Equation (2) and the isotherms of Fig. 5 can also be used to estimate devia- 


tions from ideal mixing on the individual sites from systematic deviations of the 
data from the isotherms. The excess free energy contributions for such a source 
should have the same effect proportionally as temperature variations. It is 
possible that some of the deviations from he data points from the isotherm of 
Fig. 1 are attributable to this source, but th^ uncertainties in the data preclude 
any quantitative estimate of the effect at this time. 

EQUATION OF THE EXCHANGE KINETICS 

Although no kinetic data exist for the exchange reactions in garnets, it is 
interesting to write down the phenomenological form of the rate equations. This 
may be done if we assume as did Dienes 18 that the rate controlling process is the 
exchange rather than some other process such as nucleation. Thus if Co is 
the concentration (in moles or particles per cm 3 ) of the total tetrahedrally and 
octahedrally coordinated ions we obtain the following equation 16 for the rate 
of change of concentration of one of the ions on the tetrahedral site: 


dX 


- - C K rfa X f , (1 - X a ) 


dt 5 


K arf : aH X a (1 - X ;J ) 


( 5 ) 
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This is related to rate of change on the other site as 


1 . 

2 dt fit 

K ad and K (la are the specific rate constants for exchange in the two directions. 
It is apparent that when (dX d dt)= 0 equation (5) reduces to equation (1). The 
important point to be made here is that although $ (Ja /<p ad may approach unity 
with essentially ideal mixing on the stable lattice sites, the individual 0's may 
deviate widely from unity. The reason for this is that these factors embody the 
thermodynamic properties of the transitional state and mixing in this state 
might be far from ideal. Of course, if mixing were also ideal in the transitional 
state the * 's would reduce to unity and equation (5) would assume its simplest 
form. 
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Y 3 Ga x Fe 5 _ x 0 12 



Figure 1 — Di stri but ion of Ga^ 4 and re^ among the tetrahedral and octahedral sites of 
Y Ga Fe,_ 0 1? The equilibrium constant K J is also shown. 
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